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Abstract: A one-dimensional driven lattice gas with disorder in the particle hopping prob- 
abilities is considered. It has previously been shown that in the version of the model with 
random sequential updating, a phase transition occurs from a low density inhomogeneous 
phase to a high density congested phase. Here the steady states for both parallel (fully 
synchronous) updating and ordered sequential updating are solved exactly and the phase 
transition shown to persist in both cases. For parallel dynamics and forward ordered se- 
quential dynamics the phase transition occurs at the same density but for backward ordered 
sequential dynamics it occurs at a higher density. In both cases the critical density is higher 
than that for random sequential dynamics. In all the models studied the steady state ve- 
locity is related to the fugacity of a Bose system suggesting a principle of minimisation of 
velocity. A generalisation of the dynamics where the hopping probabilities depend on the 
number of empty sites in front of the particles, is also solved exactly in the case of parallel 
updating. The models have natural interpretations as simplistic descriptions of traffic flow. 
The relation to more sophisticated traffic flow models is discussed. 
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1. Introduction 

The asymmetric simple exclusion process (ASEP) is an archetypal example of a driven 
diffusive system |^ for which analytical results are possible, particularly in one dimension 
The model comprises particles which hop stochastically in a preferred direction with hard 
core exclusion imposed. The model has a natural interpretation as a simplistic description 
of traffic flow on a one lane road and indeed forms the basis for more sophisticated traffic 
flow models 0,0. In particular one may cite variations of the model proposed originally by 
Nagel and Schreckenberg §, 0, |, |, |T0|, [11 . 



However, a basic difference between the original ASEP and traffic flow models lies in the 
updating scheme. In the mathematical literature the ASEP is usually defined in continuous 
time or, equivalently for the purposes of simulation, by a random sequential updating scheme 
where for each update a particle is selected at random. In contrast, when simulating traffic 
flow parallel updating is usually employed for reasons both idealistic — parallel dynamics 
provides a perhaps more faithful representation of real traffic — and pragmatic — parallel 
dynamics yields economy of random numbers. 

For random sequential dynamics a relative wealth of exact results on the ASEP are 
now available |0, |^, 0, 0, |T^, ^ in particular the steady states of various 

models have been constructed using a matrix product ansatz |l^ |2^, ^ |2^, |4|. This 
technique has been extended to a sublattice parallel updating scheme [[ 



H, 13] and, in 
i^. However, 



the case of open boundary conditions, to an ordered sequential scheme pq , ^ 
for fully parallel dynamics only a few exact results are known 30]. 

In the present work we determine exactly the steady state of an ASEP with disorder for 
both parallel dynamics and ordered sequential dynamics, in the case of periodic boundary 
conditions. The disorder takes the form of quenched random hopping probabilities assigned 
to each particle. Other types of disorder and inhomogeneities such as random rates asso- 
ciated with lattice sites ||3^, slow defect sites |^ or slow defect particles ||3^, ^ 



have been considered. The random sequential version of the disordered model considered 



here has previously been studied ||3g, ^ ^ and a transition shown to occur between a low 
density inhomogeneous phase, where a tailback forms behind the slowest particle, to a high 
density congested phase. 

In the present work it will be demonstrated that this transition persists under parallel and 
ordered sequential updating. The difficulty in obtaining exact results for parallel dynamics 
and ordered sequential dynamics lies in the construction of the transfer matrix. Using a 
technique inspired by ref. we explicitly construct the transfer matrices in a convenient 
form that allows the steady states to be demonstrated. 

From the point of view of traffic flow the phase structure of the disordered model is of 
interest since the disorder induces emergent jams at low densities, whereas in other traffic 
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flow models where a phase transition occurs it is in the high density phase that the jams 
emerge j^, |ll], |12|. The coarsening of the resulting jams has also been studied [^0|, |10|, ^ ^ 



From a theoretical viewpoint it was shown in the random sequential case that the transition 



has a strong analogy with Bose condensation and that the steady state velocity of a 
particle was equivalent to the fugacity of an ideal Bose gas. This analogy will be pursued 
here for the parallel and ordered sequential cases and it will be shown that the steady state 
velocity remains related to the fugacity of a Bose system. 

It is of interest to determine whether distinct updating schemes can produce different 
behaviour. It turns out that the the value of the critical density may depend on the updating 
scheme. It will be shown that the critical density is highest, implying the throughput most 
efficient, for a backward ordered updating scheme where the updating sequence is opposite 
to the direction of flow of particles. In contrast, the random sequential updating scheme 
yields the lowest critical density. 

Another key difference between traffic models (e.g. the model of |P) and the ASEP 
is that in traffic models particles may move a distance greater than one lattice spacing, 
thus implying that the dynamics is not nearest neighbour in the sense that cars are aware 
of cars several lattice sites ahead. At present analytical results are not generally available 
for the case of hopping more than one lattice site, although a step towards this goal has 
been made On the other hand dynamics of range greater than nearest neighbour may 
be mimicked by letting the probability with which a particle hops forward depend on the 
number of empty lattice sites in front of it. Indeed, the concept of a braking distance 
furnishes a natural interpretation for hopping probabilities which increase as a function of 
the empty space in front. 

In the present work we solve the steady state of a generalisation of the dynamics where 
the disordered hopping probabilities are dependent on the empty space in front, restricting 
our attention to parallel updating. 

The paper is organised as follows. In section 2 the model and updating schemes consid- 
ered are defined and exact expressions for the steady states presented. The proof of these 
expressions is somewhat technical and is deferred until section 4. Prior to that, in section 
3, the phase transition and the analogy with Bose condensation is analysed. In section 
5 the generalisation of the model to the case of space dependent hopping probabilities is 
considered and in section 6 conclusions are drawn. 
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2. Model Definitions and Steady States 

In this work we study asymmetric exclusion models where the particle hopping probabilities 
are quenched random variables. We consider M particles labelled {/x = 1, . . . , M} hopping 
on a one dimensional lattice of size N sites labelled {i = 1, . . . , N} with periodic boundary 
conditions (site N+i = site i). The random sequential version of the model was considered in 
p7| . We now define three distinct variants of the model according to the following dynamics. 

Parallel Dynamics: at each time-step all particles {/i = 1,...,M} simultaneously at- 
tempt to hop forward each with its own probability p^. A hop is only carried out if the 
target site was empty before the update. Since no backwards hops are permitted the ques- 
tion of what would happen if two particles simultaneously attempted to hop onto the same 
site does not arise. 

Ordered Sequential Dynamics: in this case a time-step corresponds to updating each 
particle in a fixed sequence. At each update the relevant particle fi attempts a hop forward 
with probability p^. We consider two sequences for the updating 

Forward updating in which the order is 1,2,... M. 
Backward updating in which the order is M, M — 1, ... 1 

Since particles cannot overtake, the sequence of particles is preserved in all three variants. 
2. 1 . Expressions for the steady states 

In order to express the steady state of each variant of the model defined above, we consider 
the weight F{ni,n2, . . . iUm) of a configuration comprising particle 1 followed by Ui holes 
(empty lattice sites); particle 2 followed by n2 holes and so on. The weights are related to 
probabilities P/v({n^}) via a normalisation Zjsi^m defined by 

PN{{n^}) = F({n4)/Z^,M. (1) 

We now present exact expressions (to be proven in later sections) for the steady state of 
the three variants. In each case the steady state weights have the factorised form 

M 

n{n,}) = n fM (2) 
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The different variants have different expressions for fi_i{n^) as follows. 
Parallel Dynamics 
In this case 



1 - p^) for 



n, 



for > 



(3) 



which may be rewritten employing the usual Heaviside function as 



Pa* 



Pa.) 



(4) 



Note that the pure case, which has P/^ = p V/i recovers a result obtained in [§] through a 

combinatorial argument. 

Forward Ordered Sequential Dynamics 

To define a steady state for ordered sequential dynamics one must specify the point in the 
update sequence to which the steady state refers. For forward updating we define the steady 
state as that after the update of particle M (the final update of the sequence). That is, the 
steady state weights are those for finding the system in a given configuration after the last 
update of a sequence, and before the first update of the next sequence. In this case 



1 - Pm 

. Pm , 



1 -Pm' 



d{n^) for /iT^M 



(5) 



Backward Ordered Sequential Dynamics 

In this case the steady state again refers to the weight for finding the system in a given 
configuration after the final update of the sequence, this time the update of particle 1. 



/m(^a/) 



;i - Pm) + 



for /i 7^ M 



1 - Pm 
Pm 



■riM 



(1 - Pm) 



0{nM) 



(6) 



Before proving the steady states (||,^|D we discuss the implications of the form (Q). 
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3. Phase transition and analogy with Bose condensation 

A steady state of the factorised form gives rise to the possibihty of a phase transition. 
In the analogy with Bose condensation was made. Let us review briefly the qualitative 
aspects of the phase transition and the analogy. 

The two phases exhibited by the model are a congested phase that exists at high density 
and an inhomogeneous phase that exists at low density. In the congested phase, the velocity 
of particles is limited by the availability of empty sites, whereas in the inhomogeneous phase 
the velocity is limited by the hopping rate of the slowest particle. Thus, the inhomogeneous 
phase may be pictured as comprising two regions: a tailback behind the slowest particle (a 
high density region) and empty space in front of the slowest particle (a low density region). 
The analogy with Bose condensation is to think of the empty sites as bosons and the state 
of a boson as determined by which car it is immediately in front of. Then in the congested 
phase the bosons are thinly spread over all the Bose states (i.e. there are uniformly small 
gaps between the particles) whereas in the inhomogeneous phase the bosons are condensed 
in front of the slowest particle (there is a large empty space in front of the slowest particle). 

Before exploring the analogy with Bose condensation we shall discuss the mathematics of 
how calculations are performed. We focus first on the case of parallel dynamics and extend 
the results to ordered sequential dynamics in section 3.3. 

To calculate quantities of interest, in particular steady state averages, one must first 
have at hand expressions for the normalisation Z^^m of the weights (||) defined through (|l|). 
Due to the constraint of — M holes the normalisation may be written 

M 

Zn,m = ^ ^Y. "m- {n-m) n (7) 

ni,n2...nM fJL=l 

This quantity may be considered as the canonical partition function of a thermodynamic 



system and in the standard way |^ it may be written using an integral representation of 
the delta function as 

^^'^ = / 1^ Z{Z) , (8) 

where the grand canonical partition function Z is given by 



M 

z{z) = n 



n„=0 



(9) 



For large N, M (||) is dominated by the saddle point of the integral and z*, the value of z 
at the saddle point, is the fugacity. 
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We now calculate the velocity f , defined as the steady state average of the rate of hopping 
of a given particle. Since the particles cannot overtake each other, the velocity is the same 
for any particle fi. Taking /i = 1 one finds 



(Af-A/+1) 



oo 

PiE(^*r/iK) 



ni=l 



M 

n 



n„=0 



ni =1 



ni=0 



for large N, M 



1 + z* 



(10) 

(11) 
(12) 



where the last equality results from performing the geometric series obtained when (^ is 
inserted in (pJ]). 

In the following the thermodynamic limit is defined by 



N ^ oo with M = pN , 



(13) 



with the density p held fixed. In order to determine the fugacity one uses the condition 

d\nZ 



N-M = z- 



dz 



(14) 



which is the saddle point condition for (|^) (or equivalently the condition that in the grand 
canonical ensemble the average number of empty sites is — M). Using (^ one finds that 
(U) becomes 

;i -P/^)^ 



M 

z{z) = n 

11=1 L 



1-Pm + 



- (1 -Pm)^. 



and (0) yields 



N-M 



M 

E 



1 + - (1 -Pm)^ 



(15) 



(16) 



Now using the relation between z* and v ([T^), one obtains in the thermodynamic limit 

1 *^ 1 

1-P = ^(1-^)t7E Z -. (17) 



3. 1. The disorder average 

In the following we shall consider the particle hopping probabilities as quenched random 
variables drawn from a common distribution |38| 



:i8) 
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with support on the interval [c, 1). Other distributions may also be considered |3^, but the 
qualitative behaviour is determined by the power 7. 

For the disordered system one wishes to obtain properties given by a typical realisation 
of the disorder (here the particle hopping probabilities), for example the typical velocity. 
Usually in the theory of disordered systems one expects quantities such as v to self-average 
but quantities exponentially large in the system size, such as Z^^m, not to. That is, one 
expects as N ^ 00, v ^ v with probability one, where the bar indicates an average over 
the particle hopping probabilities. When working in the grand canonical ensemble one has 



directly an equation for v in the thermodynamic limit, (17). Thus, with the assumption 
that V is self-averaging, one simply has to average ([T7|) to obtain an equation for the typical 
velocity. 

Now, in the thermodynamic limit (0), the fraction of particles with hopping probabilities 
between p and p + dp converges to V{p)dp. Equally the velocity of the slowest particle 
converges to c. For the sake of clarity, however, it is convenient to assume that the slowest 
particle, taken to be particle 1, has velocity exactly equal to c. Therefore (0) may be 
replaced by 

l-p = p{l-v)I{v)+^^ (19) 

where ^ 

I{v) = I dpV{p)-^ (20) 

Jc p — V 



Since the rhs of (|1^) is an increasing function of z we deduce that the rhs of {\17\) is an 
increasing function of v. However, z was introduced in (|]) as a contour integration variable 
therefore the saddle point z* must lie between and any pole in the integrand of (H). Owing 
to this, the maximum value z can take is c/(l — c) and the maximum value of v is c. 

If the integral (^) diverges as f — > c then (|19D can always be satisfied with v < c and 

< rii > /N zero. However, it turns out that for the distribution (|I8|) with 7 > 0, I{v) is 
always finite as v c from below. Therefore when 

which holds for p below a critical value pc, (|T^ can only be satisfied with u — »• c and 

< ni > /N non-zero. As (|19D expresses the constraint in the number of holes, we see that 
in this case a finite fraction of the holes must reside in front of the slowest particle. Thus, 
we have the inhomogeneous low density phase. On the other hand for p > pc, (|19D may be 
satisfied for v < c and < rii > /N zero, in which case we have the high density congested 
phase. 

By this point the qualitative analogy with Bose condensation should be apparent. How- 
ever, we first remark that the exact mapping onto an ideal Bose gas found for random 
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sequential dynamics [T^, ^] no longer holds. In the mapping is viewed as the occupation 
number of the /i**^ Bose state; applying this to Eq. H one sees that in the parallel update 
system an unoccupied Bose state is penahsed. This effective repulsion between particles 
already occurs in the pure case Secondly, for random sequential dynamics the velocity 
was exactly equivalent to the fugacity of an ideal Bose gas |^ , whereas in the present case 
although V is still an increasing function of z*, the relation is modified to (0). Nevertheless, 
the analogy with Bose condensation remains useful. The slowest particle corresponds to the 
Bose ground state and the distribution of particle hopping probabihties V{p) corresponds 
to the density of states of the Bose gas. Thus the transition to the inhomogeneous phase 
corresponds to a condensation of a finite fraction of bosons into the Bose ground state (i.e. 
a finite fraction of the holes reside in front of the slowest particle). 



3. 2. Critical Behaviour 



e = In 



(22) 



In order to analyse the phase transition we require at our disposal the asymptotic behaviour 
of (|T^ as f — c. The expansion in e where 

"c1 c — V 

- ~ as f — >• c 

_f J c 

is carried out in appendix A. The result is 

(I _|_ / C \ ^+'^ 

For 7<1 liv) =i^^ + (l + 7)(^^j r(l + 7)r(-7)e^ + 

2c 



For 7 = 1 I{v) =^_^ + 2(^) eln(e) + . 

(1 + 7)c (1 + 7)c " 



For 7 > 1 I{v) 



1 + 



7 



(1-c) (7-1) 



e + . 



(23) 
(24) 
(25) 



7(1 - c) 7(1 - c) 

we see from (p3|-p5|) that /(c) = (l+7)c/ (7(1— c)) which implies, together with the condition 
for the transition (^1]) , that the critical value of p is 



Pc = • (26) 

7 + c + 7c 

One should note that for that 7 = (a fiat distribution of hopping probabilities) the 
transition to a congested phase occurs at zero density. 

As p is increased above the critical value the velocity decreases from c according to 

' ^\p - Pc)'/'' + . . . (27) 



For 7 < 1 V ~ c — 
For 7 = 1 V ~ c — 
For 7 > 1 V ~ c — 



:i-c)^ 



+ . . . 



p2cr(i + 7) |r(-7) 

(1 - C) (p - pe) 

2p2 c I ln(p - Pc 

(1 _ c)7(7 - 1) 

p2 (1+7) [c+(7-l)(l-c)] 



(P - Pc) + 



(28) 
(29) 
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One can also consider the current or throughput J defined by 



J = p V . 



(30) 



In the inhomogeneous phase {p < pc) the current increases hnearly with p and J = cp. 
Using the expansions (|27|-p9D the behaviour of the current in the high density congested 
phase can be analysed. Just above the transition the current increases with p if 7 < 7* but 
decreases if 7 > 7* where 



7 



c + c^ 



2(1 -c) 



1 + J1+4(1-C)CV(1-C + C2)2 



(31) 



Thus the maximum current is achieved in the high density phase if 7 < 7* and exactly at 
the phase transition if 7 > 7*. However it does not appear that 7* is universal because it 
depends on c i.e. it depends on the particular choice of V{p). 



3. 3. Ordered Sequential Dynamics 

For forward sequential or backward sequential dynamics the velocity may be defined as 
the probability that the first particle in the updating sequence hops forward. It is easy to 
check that (|T2p is still the correct expression for the velocity in both these cases so that the 
velocity is always related to the fugacity in a simple way. By performing the sums in (|^) 
one obtains for forward updating 



A/-1 

n 



Ziz) 

and for backward updating 
Z{z) 



1 



Pm 



Pm 



AI-l 

n 

11=1 



1 



Pm + 



:i 



Pm)z 



Pm)z 



(32) 



Pm - (1 -Pm)z 



(33) 



For forward updating the saddle point value of the fugacity z and hence the velocity 
will be the same as as for parallel updating, implying the same critical density. However, 
for backward updating (0) gives a different saddle point value of z to (plSf ), and it can be 
checked that (P^Ql) is modified to 



l-p 



(1 — v) I{v) — V 



(34) 



givmg 



7 

c + 7 



(35) 
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We see that for backward updating pc is greater than for parallel and forward updating 
which share the same pc given by (0). Since the velocity in the low density phase is greater 
than that in the high density phase, this means that backward updating yields the best 
throughput. One might expect this because backward updating increases the chance that the 
site in front of the next particle to be updated has just been vacated. One can also calculate 
the critical density for random sequential dynamics [0, ^ to find pc = 7(1 — c)/(c + 7) 
which is lower than both ( pSj) and (|26|). This implies that random sequential updating gives 
the poorest throughput. 



4. Proof of Steady State for Ordered Sequential and Parallel Dynamics 

In all three variants of the model considered here the updating rules comprises a determined 
time-step. Therefore a transfer matrix may be used to express the condition for the steady 
state weights as 

Y^T{CX')F{C) = F{C) , (36) 
c 

where T(C, C) are the components of the transfer matrix. In (|36D T(C, C) is the proba- 
bility of going from configuration C to C in one time-step and r(C,C) is the probability 
of remaining in configuration C after a time-step. A configuration is specified by the hole 



occupation numbers C = {rii, . . . , njv/}- We see from ( |36| ) that the steady state weights form 
an eigenvector of the transfer matrix with eigenvalue one. 

In the construction of the transfer matrix that follows, it is convenient to use an oper- 
ator notation. To each configuration {rii, . . . ^hm} is associated a vector | rii, . . . ,nM), to 
be referred to as the configuration vector. The configuration vectors | rii, . . . ^um) form an 
orthonormal basis for the vector space of configurations. We write the weight of a configu- 
ration as 

F(ni,...,n^) = (F| ni,...,n^) (37) 

Thus, the the weight of a configuration {ni,...,nM} is the component in the direction 
(ni, . . . , of the bra vector {F\. 

With this notation we may rewrite (PB|) as 



(F|f|ni,...,n^) = (F|ni,...,n^) (38) 

where T is an operator acting on the space of configuration vectors defined above: T acting 
on a given configuration vector generates the possible configuration vectors before an update 
multiplied by the appropriate transition probabilities. We shall refer to T as the transfer 
matrix as well as its components T(C,C'). Our task now is to construct T. 
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4. 1. Construction of the Transfer Matrices 

In this subsection we construct the transfer matrix for all three variants of the dynamics. 
The desired form is of a trace of a product of matrices, each of which contains as elements 
operators acting at the relevant site |3^. The technique is most directly illustrated in the 
case of ordered sequential dynamics, which we consider first in this subsection. 

The transfer matrix for a full time-step of ordered sequential dynamics may be written 
as an ordered product of M operators h^_i ^ corresponding to the update of each particle ix 
in sequence. Recalling that T acting on | ^i, . . . ,n^) generates the possible configurations 
leading to | ni, . . . , n^), one deduces that to correctly generate these configurations we first 
act on I ni, . . . , n^) with an operator corresponding to the last update of the sequence, then 
with an operator corresponding to the second last update and so on, so that for the forward 
updating 

M 

= n K-i A* • (39) 
By similar reasoning, for backward updating the transfer matrix Tb is 

M 

Tb=W hM-f, M-fi+i ■ (40) 

In both cases the operator ^ is given by 

h^_i ^ = 1 - Pf^a^al + Pf,al_^a^ (41) 

where is a raising operator and aj^ is a lowering operator acting on the vector space 
spanned by |ni, . . . , um >■ 

a^|ni, . . . ,n^, . . . ,nM > = \ni, . . . + 1, . . . ,nM > (42) 
al\ni, . . . ,n^, . . . ,nM > = \ni, . . . - 1, . . . ,nM > if rif, > (43) 

= if n^ = 0, (44) 

and one may verify that 

ii\ni, . . .,nM >= {1-Pf,9{n^))\ni, ...,nM> +Pf,0{n^^i)\ni, . . . ,n^^i-l,n^+l, ...,nM> 

correctly giving the configurations leading to \ni, . . . ,nM > after the update of particle fi. 

To proceed further we employ techniques used in the study of integrable models . 
First it is easy to check that (^) may be rewritten in two ways as 

h = (1. f I-WM (45) 



(1 -p^a^at, p^ttf,) ( 1 



(46) 
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Equations (^Sf ) and (|4q) are to be read as a scalar product in an auxiliary space (see below). 
On inserting (|45|) into (^), Tp becomes 



% 4/) 



n 



1 - PMauO^M 

Pmo-m 



(47) 



which may be rewritten as 
Tp = Trace 




1 - p^,af,al al - p^^a^a^a^^ \ ( PMauali (1 - Pm)4^ 



Pt^o-t^ 



Pm^m 



Pm'^ 



(48) 



To obtain (48) from (^) it may be checked that for arbitrary commuting operators Xj and 
yi and arbitrary operators Zi, one has the identity 



{xi, X2] 



yi y2 




Trace 



yi y2 
ys 2/4 



X\Z\ X2Z1 
XIZ2 X2Z2 



(49) 



We stress here that (^8]) is merely a convenient way of writing the sums of products of 
the operators 1, a^, aj^ which form the transfer matrix. The two by two matrices appearing 
in (^8]), whose elements are made up of the operators 1, a^, aj^ may be thought of as acting 
in some auxiliary space. The trace is carried out in this auxiliary space and not in the space 
in which the operators 1, a^, aj^ act. 

In a similar fashion one can construct the transfer matrix for backward updating, this 
time using (^61) in (EOf). One obtains 



Tb = Trace 



1 



Pmo,m<i-m Pmo-m 




PM-hO,M-h 
PM-u^ 



(50) 

Finally, to construct the transfer matrix for parallel dynamics Tp one should first realise 
that it is closely related to Tp since for forward updating each particle is unaffected by the 
results of previous updates in the sequence, except for the final particle M. After a little 
reflection it can be confirmed that 



fp = Trace 



M 

n 



1 - p^a^al 

P/iClfj, 



,t _ 



p^a^alal 



p^a^a 



(51) 



4. 2. Proof of the Steady State 

It was stated in section 2.1 that the steady state is given by (H), with f^iji^) obeying (^) 
for parallel dynamics, (^) for forward sequential dynamics and (P) for backward sequential 
dynamics. The proof is quite similar in all three cases. We shall describe it in detail first 
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for the case of parallel dynamics where the transfer matrix has the form (plf ). At the end 
of the subsection we shall return to forward sequential then backward sequential dynamics. 

We first note that the vector space spanned by \ni, . . . ,n^) is a tensor product of M 
spaces each with basis vectors \n^ > where = ■ ■ ■ oo 



(52) 



Thus the action of the transfer matrix for parallel dynamics (^) on ( |52D may be written as 



Tp\ni, . . . , um) = Trace 



AI 

n 



(53) 



where the product, in fact, indicates a tensor product in the configuration space (as well as a 
usual product in the auxiliary space of two by two matrices) and the operators 1, a^, a|j act 
on the space spanned by \n^). To keep the notation light we have not explicitly indicated 
these matters, but the meaning is clear. 

We are now in a position to prove the expression for the steady state given by (0) and 
(I). If the steady state is of the factorised form (Q) we have 



(^1 = E [ ] ® [ /2K)(ri2| ] ■ ■ ■ ® [ fMMiuMl 



(54) 



and one finds from ( [53|j5^) that 

{F\Tp\ni, . . . , Um) = Trace 

where 

B,{0) - 
B,{1) - 



M 

n 



n B,{n,) 



(55) 



and for n> 1 Bi^i{n) 



( fM o\ 

/ (l-p,)/,(l) /,(0) \ 
\ P^,f^,{n + l) Pt,ff,{n) 



(56) 
(57) 
(58) 



Inserting the expressions for /^(n) given in (^, and employing the Heaviside function we 
may rewrite (^^-|58D as 

B,{n) = Un)A^{n) (59) 

where 

/ 1 — n Mini n H(n\ \ 

(60) 



1 - p^9{n) p^9{n) 
l-p^9{n) p^9{n) 
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To prove the steady state (0) we must show that {F\Tp\ni, . 
on inserting (p9D into (|55|) , we see that it remains to show 



Trace 



M 



riAd) = Ua=lf^^M■ Thus, 



(61) 



In order to do this we seek a similarity transformation 



A„(n) = L,,A,,(n)R„j.'i with R„L„ = 1 



that puts A^{n) into the form 



1 X 





(62) 



(63) 



where x is any number. The trace of a product of such matrices ( |63D is clearly unity thus 
satisfying (|6T|) . A similarity transformation which fulfills (|62| , |63D is straightforward to 
construct. Taking 



1 P^^ 



it is easy to verify that R^L^ = 1 and 



(64) 



A^{n) = Lf,Af,{n)Rf,+i 



1 Pt^+i - Pf^0{n) 




(65) 



hence @) is proven. 

In a similar manner the steady state for forward updating may be proven. Here we only 
provide a few key points. Following the steps (|53|) through to (^O]), now using (|48|) as the 
transfer matrix one finds that ( |60D holds for fi < M and 



1 - PMO{n) pM0{n) 
1 - Pm Pm 



(66) 



The similarity transformation (|6^) then gives 



-4m (n) 



1 Pi - - pm{1 - PM)d{n) 







pm(i - e{n)) 



(67) 



whilst the other Aij_{n) have the form (p5|). The trace of a product of these A/x('^) again 
gives unity. 

For backward updating, using (|50D for the transfer matrix, one obtains for /x < M 



1 - p^,d{n) l-Pf, 
Pf^O{n) Pf, 



(68) 
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and for fi = M 
A suitable similarity transform is now 



AM{n)=('-P-f:^^ ^-P^f-)] (69) 
V PMO{n) PMO{n) I 



4K) = Rl+iM^,)Ll (70) 



where i?^, are still given by (|64|). One then obtains for /i < M 

'^'^''^ ^ ( -pI- - PMn) - d{n)) ) ^^^^ 

and 

^M{n) =( ^ n V (72) 

It is easy to convince oneself that the trace of a product of matrices of form (^) with a 
single matrix of form ( [72D yields unity, thus completing the proof of the steady state for 
backward updating. 

5. Generalisation to hopping probabilities dependent on empty space in front 
of particles 

In this section we consider hopping probabilities that depend on the number of empty sites 
in front of a particle. For this dynamics with random sequential updating, it is known 
that the steady state is given by a product measure Here for the disordered case we 



restrict our attention to parallel updating although it is straightforward to generalise to 
ordered sequential updating. Thus for parallel updating, at each time-step particle fi will 
hop forward with probability p^(n^), where we recall is the number of empty sites in 
front of particle /i so that clearly p^(0) = 0. In a traffic model suitably chosen hopping 
probabilities could mimic the effect of deceleration when another car is encountered. 

It turns out that the steady state again has the form (|^). It will be shown in this section 
that 

f,{n) = (l-p^(l)) for n = 
^-pMk=i PA^) 

In order to prove (|73|) we generalise the operators used in (^T|-^) by defining new operators 
a^{n),b^{n): 

a^{n)\ni, . . . ,nf„ . . . ,nM > = |ni, . . . , + 1, . . . , um > if n,, = n (74) 

= otherwise (75) 

bf,{n)\ni, . . . ,nf„ . . . ,nM > = . . . , ra^, . . . , riAf > if n^ = n (76) 

= otherwise . (77) 
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The corresponding operator to (|4lD is now given by 

oo oo 



n] 



(78) 



n=l 



n=l 



where aj,_]^ is still defined by (^^4|). Using the same procedure outlined in section 4 the 
transfer matrix (for parallel dynamics) may be constructed as 



tp = Trace 



M 

n 



V EZiP.in + 1) a,in) [j:n=iP,in + l)a^(n)] 4 



(79) 



Using this transfer matrix and assuming that the steady state is of the form (^, in analogy 
with section 4 one arrives at (E^) where now 



^(n)(l - p,{n)) f,{n - 1)(1 - p^{n - l))e{n) 
f^{n + l)p^(n + 1) f^{n)p^{n) 



(80) 



and it should be borne in mind that p^(0) = 0. Inserting (|73|) one finds that it remains to 
show condition ( pT] ) is satisfied where this time 



A, 



I 1 -p^(n) p^{n) 



(81) 



It is easy to check that inserting _R^(?2^)L^(n^) (= 1) before each y4^(ri^) in the product 
(|6l|), where 



n„ 



V 1 -(1 -PmK)) , 
yields the required similarity transformation. 



and Lfj,{n) 



l-Pi,{nf,) P^in,,) 
1 -1 



(82) 



6. Conclusion 

In this paper it has been shown that analytical results may be obtained for the asymmetric 
exclusion model with parallel dynamics which forms the basis for many discrete models of 
traffic fiow. Ordered sequential dynamics have also been treated using similar techniques. 
Previously, the steady state for the situation where at each parallel time-step each particle 
attempts forward with probability p had been solved Here we have generalised that 
model to one where at each time-step each particle fi attempts a hop forward with probability 
p^(n^), an arbitrary function of the number of empty sites immediately ahead of particle 
fi. The method used to prove the steady state ( |T5D was to construct the transfer matrix 
using a technique inspired by p9 . 
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A detailed analysis was made of the case where did not depend on n^, and it was 
shown that a phase transition analogous to Bose condensation occurs, as it did for a random 
sequential updating [^]. In particular it is interesting to note equation ([I^ ) which shows 
the velocity in the particle system is intimately related to the fugacity (or equivalently 
chemical potential) when the system is viewed as a Bose system. Since in (single component) 
equilibrium systems the thermodynamic phase is determined by minimising the chemical 
potential, the suggestion is that in simple traffic flow models a principle of minimisation of 
velocity pertains. 

A simple illustration of this principle is the example of a single slow particle with hopping 
rate c while all the other particles have hopping rate p ||3^. Then in a phase where there is 



a tailback behind the slow particle and empty space in front, the velocity is c whereas in a 
congested phase the velocity should be given by ( 1 — y^l — 4pp(l — p) )/(2p) Choosing 
the phase with the minimum velocity yields the critical density Pc = {p — c)/{p — c^). This 
expression is, in fact, exact as can be checked by using the results of section 3. It is important 
to ascertain whether a principle of minimisation of velocity holds in more complicated traffic 
flow models. 

It was also interesting to note that parallel dynamics and forward sequential dynamics 
have the same thermodynamic behaviour e.g. the same critical density. One expects this 
since a time-step of forward updating and a parallel time-step only differ in last update 
of the sequence in forward dynamics. However backward updating has a distinct critical 
density as does random sequential updating. Comparing the three critical densities reveals 
that the highest is for backward dynamics, the lowest is for random sequential dynamics 
and parallel dynamics lies in between. This conflrms that backward dynamics is the most 
efficient updating scheme in terms of throughput and random sequential is the worst. 

It would be interesting to explore other realisations for the case of general p^{np) con- 
sidered in section 5. In particular it should be possible to analyse the effect of a braking 
distance for each particle. 

The present work narrows the gap between bona fide models of traffic flow, such as that of 
Nagel and Schreckenberg or more sophisticated models 0, and the simple particle hopping 
models for which exact results are possible. In order to close the gap further, analytical 
results are desirable for models where particles can hop more than one lattice site at each 
update. Some progress has already been made towards this . 



It would also be of interest to analyse more complicated properties of the system such 



as relaxation to the steady state, as studied numerically in |3^, |41[], or else the diffusion 
constants of particles O, 111 
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Appendix A: Evaluation of the Integral ( P0|) 

In this appendix we derive the expansion (p3|-|25|) of the integral (pO]) given by 



By first defining 



1|) becomes 



Uv) = — —— / dp — — . (Al) 



V = cexp(-e) , (A2) 



Tf \ 7 + 1 f^-"", cexp(-e)p^ 



and we may follow an analysis similar to that of Robinson We make the Mellin 

transformation 

I{s) = / de I{e)e^-' , (A4) 



n+l -p/ N 

1 r S 



p + cj {n + iy 



(A5) 



which may be expanded to yield 
where r(s) is the usual Gamma function defined by 

roo 

T{s) = / de exp(-e) e*"^ . (A6) 

The inverse transformation is 

rx+ioo 

/ ds I{s)e-' , (A7) 



2 rX+tOO 

2Td 



where x is a real constant chosen so that the contour of integration is to the right of any 
pole. The integral can be evaluated by closing the contour and using the calculus of residues. 
The analytic structure of 

X{s) = g,{s)T{s) , (A8) 

where 

g Js) = , / dpp^y^i — , (A9) 

is as follows: T{s) has simple poles at s = —n where n = 0, 1, . . . oo with residues (— l)"/n! 
and g-yi^s) has a simple pole at s = —7. To evaluate the residue of ^'7(5) at s = —7 we note 
the small p behaviour of the sum involved in ( |A9| ) 

00 / ^ 
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Thus, one obtains 



and 



/(e) = (1 + 7) (t^) r(l + 7)r(-7)e^ + £ ^ 

'^^ n=0 



The expansion (p3| , p5D is the first two terms in (|A12|) where it can be computed that 

1 + 7 c 



7 1 — c 
1 + 7 c 
7 1 — c 



1 + 



7 c 



7 — 1 1 — c 



(All) 
(A12) 

(A13) 
(AM) 



For 7 — s> 1, the singularity in ( [A14| ) cancels with the singularity in the first term in the rhs 
of ( |A12| ) and one obtains 



/(e) = g,{0)+^m (1 + 7) (^_- j r(l + 7)r(-7)e^ - ^7,(-l)e 



+ 



9^(0) + 2 [^^J elne + 0(e). 



(A15) 
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